Klein Backscattering and Fabry-Perot Interference in Graphene Heterojunctions 
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We present a theory of quantum- coherent transport through a lateral p-n-p structure in graphene, 
which fully accounts for the interference of forward and backward scattering on the p-n interfaces. 
The backreflection amplitude changes sign at zero incidence angle because of the Klein phenomenon, 
adding a phase tt to the interference fringes. The contributions of the two p-n interfaces to the phase 
of the interference cancel with each other at zero magnetic field, but become imbalanced at a finite 
field. The resulting half-period shift in the Fabry-Perot fringe pattern, induced by a relatively weak 
magnetic field, can provide a clear signature of Klein scattering in graphene. This effect is shown 
to be robust in the presence of spatially inhomogeneous potential of moderate strength. 
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The electron system in graphene features surprising 
connections with the relativistic quantum mechanics of a 
massless Dirac particle in external electric and magnetic 
fields []], 0, & 0- I n particular, charge flow in lat- 
eral p-njunctions in graphene is described by Klein scat- 
tering 0, 0, 0] , exhibiting perfect transmission through 
the barrier at normal incidence, regardless of the barrier 
characteristics, and a barrier-dependent, finite reflection 
coefficient at non- normal incidence Q. On the exper- 
imental side, while transport properties of the first p- 
n junctions fabricated in graphene were dominated by 
disorder @, 0, [To[ , two recent papers Jj], [l2[ report on 
observations of a contribution to the conductance consis- 
tent with expectations for ballistic transmission through 
p-n interfaces. To extract this "intrinsic" contribution, 
which is relatively small in magnitude, one must account 
for the screening of the gate potential [13[ and the effects 
of disorder 

What other features, besides collimated transmission, 
may serve as an experimental signature of Klein scatter- 
ing? Here we focus on the characteristic behavior of the 
reflection amplitude, which exhibits a jump in phase by 
7r when the incidence angle a is varied from positive to 
negative values (see FigJT]). The sign change occurs at 
normal incidence because the reflection amplitude van- 
ishes at a = 0. Below we show that this phase shift, 
which is fundamental to Klein scattering, could serve as 
a hallmark of Klein physics in graphene. 

The backreflection phase can be detected from inter- 
ference of electron waves scattered on two parallel p-n 
boundaries in a p-n-p structure. Transmission in this 
system, described by the Fabry-Perot (FP) model, ex- 
hibits periodic dependence on the phase AO gained by an 
electron bouncing between the p-n interfaces (see Fig IT]), 



A(9 = 2(9wkb + A6>i+A6> 2 , 



(1) 



where #wkb = \ Ji p x {x')dx' is the WKB phase and 
A#!(2) are the backreflection phases for the interfaces 1 
and 2, exhibiting a 7r-jump at zero incidence angle a. 
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FIG. 1: Schematic of electron transmission through a p-n- 
p structure at zero (a) and finite (b) magnetic field B. The 
amplitude of backeflection changes sign at the incidence angle 
a = 0. At finite B, for the trajectories satisfying condition J2} 
the angles of incidence on both p-n interfaces are of the same 
sign. This adds a phase ty the electron phase accumulated 
between reflections, resulting in half a period shift of Fabry- 
Perot interference fringes, (c) Weak spatial inhomogeneity 
does not alter the relative sign of incidence angles. 

As illustrated in FigHJ the contribution AOi + AO2 to 
the net phase can be altered by a magnetic field. At zero 
B the incidence angles at interfaces 1 and 2 have oppo- 
site signs, and thus the jumps in AOk^) cancel. However, 
for curved electron trajectories at a finite £>, the signs of 
the incidence angles can be made equal. Indeed, because 
of translational invariance along the p-n interface, in the 
presence of a magnetic field the ^-component of electron 
kinetic momentum varies in space as p y (x) = p y — eBx, 
where p y is the conserved canonical momentum compo- 
nent that labels different trajectories. For the incidence 
angles at interfaces 1 and 2 to be of equal sign, p y {x\) and 
P y {x2) must have opposite signs, which happens when 



eBL/2 <p y < eBL/2. 



(2) 



In this case the net backreflection phase AOi + AO2 in 
(PQ) equals tt. As we shall see, the backreflection phase 
manifests itself as half a period shift of the FP fringe 
contrast. This phase shift develops for the field strength 
such that the range ([2]) exceeds the Klein collimation 
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range in which the p-n interface is transparent, eB > 
Ap/L (a similar condition appeared in the discussion of 
magnetoresistance in Ref.|4[). 

One useful characteristic of the phase jump in back- 
reflection is that it is less momentum-selective than col- 
limated transmission. A potential difficulty, however, is 
that the interference of scattering on two p-n interfaces 
can be sensitive to disorder. Below we will investigate 
the dependence of the FP contrast on magnetic field in 
the presence of large-scale spatial fluctuations. We find 
that, while the FP fringe contrast is suppressed, the 1/2- 
period shift, controled by backreflection, remains surpris- 
ingly robust. Even at a relatively high disorder strength, 
when the FP contrast is strongly reduced and the fringes 
become aperiodic, the 1/2-period shift induced by the 
magnetic field remains clearly discernable. 

Here we shall model the gate potential by a parabola 
U(x) = ax 2 — £, which at e > creates p-n interfaces at 

x = d=x e , x £ = \fefa. (3) 

The energy e is controled by a top gate, with e = aVt g 
(in Ref.pjJ, a ~ — ^ ). The curvature parameter a is de- 
termined by the width of the top gate L tg and its height 
h above the graphene plane. The actual potential profile 
may be nonparabolic (see Ref. for modelling of screen- 
ing effects), however this difference should not matter for 
FP interference, occurring when the separation between 
p-n interfaces L = 2x £ spans only few de Broglie wave- 
lengths. This is the case near the threshold in Vt g at 
which the p-n-p structure forms, where L is small com- 
pared to h and L tg . The parabolic U(x) may provide a 
reasonable approximation also for the devices described 
in Ref. j 121], especially those with narrow top gates. 

The Hamiltonian for a two-component Dirac wavefunc- 
tion, in the presence of the potential U(x), is 

H = v F a 3 p x + v F a 2 (p y - eBx) + U (x) (4) 

where 02,3 are Pauli matrices, vf ~ 1.1 ■ 10 6 m/s is the 
Fermi velocity in graphene, and B is the magnetic field. 
Hereafter we set ft = vf = 1, with the units for energy, 
length and magnetic field as follows: 

£* = (avph 2 ) 1 / 3 « 14meV, x* = Tivp/s* ~ 53 nm, 
B, = $o/27rxl « 0.24 T, ($ = h/e), (5) 

where we used the curvature a of 5eV//im 2 , obtained 
from the model potential of Ref. [ll[ fitted to a parabola. 
Because of the weak a 1//3 dependence, the estimates (|5]) 
should also apply, at least roughly, to other systems. 

To find the transmission and reflection coefficients, we 
factorize the wavefunction as i/j(x,y) = e %VvV ^(x\ and 
solve the one-dimensional Schrodinger equation 

id x tj) = (U (x)a 3 - i(p y - eBx)^) ^, (6) 

where without loss of generality we set the Fermi energy 
equal to zero. Transmission, evaluated from numerical 
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FIG. 2: Transmission of the p-n-p structure (a) at zero mag- 
netic field and (b) at a finite field of B — 0.2$o/x*. Electron 
momentum p y is measured in units of p* = vf/s*- Note that 
at finite B the fringe contrast is reversed across the parabola 
((7j) (black line), with the maxima and minima of the fringe 
pattern interchanging. This behavior is in agreement with the 
Fabry-Perot model ©, with the reflection amplitude vanish- 
ing on the parabola, and changing its sign. 

solution of Eq.(|6j), exhibits resonances as a function of 
energy and momentum, shown in Fig|2j We note a dras- 
tic difference between the results at zero 5, which are 
identical to those of Ref.[15[, and the results at finite B. 

To understand the behavior of transmission, it is in- 
structive to view the differential equation ([6]) as a fic- 
titious time-dependent Schrodinger evolution of a two- 
level system, the coordinate x playing the role of time. 
In this analogy, the system is driven through an avoided 
level crossing with splitting —2i(p y — eBx) at the crossing 
times determined by degeneracy of the diabatic states, 
U(x) = ax 2 — e = 0. This condition yields x = ±x £ , 
Eq.(j3]). A Landau-Zener transition at the first level 
crossing creates a coherent superposition of the dia- 
batic states that can interfere at the second crossing. 
This so-called Stiickelberg inteference (see [17[ and ref- 
erences therein), which can be constructive or destruc- 
tive, is described by an oscillatory function of the phase 
AO = — 2 j^_ £ U(x)dx = \ex £ gained between the cross- 
ings. The locations of interference fringes, determined 
from the condition AO = 2ixn with n = 1,2..., are 
e n = (37rn/2) 2 / 3 £*, which agrees with the positions of 
the fringes seen in Figl2l 
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For 5 = 0, the suppression of oscillations near zero 
p y can be linked to the absence of Landau-Zener transi- 
tions at vanishing level-splitting. In the scattering pic- 
ture, this is nothing else than the Klein phenomenon of 
perfect transmission at normal incidence. At finite £>, 
the oscillations disappear when one of the level splittings 
p y ±eBx £ vanishes, suppressing one of the Landau-Zener 
transitions. In terms of electron motion, this condition 
is equivalent to the requirement of normal incidence on 
either of the interfaces (|3]), giving 



Py 



(7) 



which is the black parabola drawn in Fig|2)3. Indeed, 
fringes disappear on this line; upon crossing the line, the 
maxima and minima of fringes interchange, indicating a 
7r phase shift in the phase of fringe contrast. 

A more refined description can be obtained from a qua- 
siclassical solution [l5| with position-dependent momen- 
tum p x {x) 



IU 2 (x) —py(x). The turning points, de- 
fined by p x = 0, are arranged as xy < x\ < X2 < xy, 
with X2 and xy equal 



e + 6 2 - Py + b, J e + b 2 + p y - b 



(8) 



where b = \eB and Xi^{p y ) — —x 2 (2')(—Py)- Hereafter 
we set a = 1, restoring physical units later. Remarkably, 
the conditions x\ = xy and x 2 = x 2 > , which correspond 
to one of the p-n interfaces becoming transparent because 
of the Klein phenomenon, yield a relation between p y and 
e which is identical to Eq.Q found above. 

The classically forbidden regions xy < x < x\ and 
X2 < x < X2> , where p x is imaginary, correspond to the 
Klein barriers at the interfaces 1 and 2. Denoting the cor- 
responding transmission coefficients as t\ and £2, we can 
write the net transmission of the entire p-n-p structure 
in a general Fabry-Perot form 



T(s,Py) 



ht 2 



iA6\ 



(9) 



where 7*1(2) = 1 — ^1(2) are the reflection coefficients, and 
the phase A6 is a sum of the WKB part and the phases 
of the reflection amplitudes, Eq.(pQ). 

The transmissions ~t\(2) can be evaluated in the WKB 
tunneling approximation: 



-2Im 'p x {x')dx' ^ e -X(p y -eBx £ ) 2 ^ . 



2ax £ 



(10) 



with the integral computed by linearizing U (x) near x = 
x £ . Similarly, linearizing U(x) near x = — x e , we find 
. Thus, the reflection amplitudes are 



l 2 — e -\{p y +eBx 



sgn(p y ± eBx £ ) e v 



vT 



-\(p y ±eBx £ ) 2 



(11) 



where we factored the sign, responsible for the phase 
jump, and a regular part of the phase e 26>reg , as follows 
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FIG. 3: Fringes in resistance (^3]) at B = 0, plotted in the 
units of i?* = (x*/VK)/i/4e 2 (blue line). Inset illustrates a 
n 2 / 3 scaling for the maxima and minima of R, which is con- 
sistent with the £ 3 / 2 dependence of the WKB phase (|12[) . 
Averaging over smooth potential fluctuations, described by 
a sum of a few harmonics, suppresses fringe contrast (red 
dashed line). Here we use Eq. (|14[) with \a m \ = 3s*. 

from analyticity in p y . Because the WKB treatment is 
exact for linear potentials [18[ , and the transmission ([TO]) 
is exponentially small unless \p y ± eBx £ \ ^ A -1 / 2 , the 
linearization of U (x) used to evaluate the integral in ([TP]) 
gives accurate results for the energies of interest, e > e*. 

The dispersion of the resonances in Fig|2] can be un- 
derstood from the momentum dependence of the quasi- 
classical WKB phase (for simplicity, we set B = 0): 



#WKB 



p x (x f )dx f 



where an expansion in the parameter \p y /e\ <C 1 is legit- 
imate because Klein collimation restricts transmission to 
Ap y ~ A -1 / 2 . The quantization condition 612 = n(n+ 1) 
gives the resonance energies £ n (Py) dispersing as in Fig 01 
To summarize, the FP model (|9j) is in full agreement 
with our numerical results. In particular, it explains the 
striking difference between the behavior at zero and fi- 
nite 5, as well as the phase shift of the fringe pattern, 
resulting from a sign change of the reflection amplitudes, 

Eq.dDJ. 

These results can now be applied to analyze conduc- 
tance and resistance, given by 



(13) 



where W is the width of the p-n structure (see Fig[T]). As 
illustrated in Figj3l the resistance exhibits fringes which 
obey the n 2 / 3 scaling, as expected from the phase depen- 
dence on the energy, Eq. ([T2|) . Somewhat surprisingly, 
the integral over p y in ([T3]) yields a fairly high fringe 
contrast in R. This is in agreement with the dispersion 
of resonances, Eq. ([T2]) . being weaker than the dispersion 
of transmission 1 1^2) (Py)- 

In the presence of magnetic field, alongside with the 
overall increase in resistance, we observe that the fringes 
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FIG. 4: Resistance ([13)) of the p-n-p structure as a func- 
tion of magnetic field and energy. The quantity plotted is 
log with i?* = (x*/VK)/i/4e 2 . Note that the resistance 

minima at high B, marked by dashed lines, are shifted by half 
a period relative to those at zero B. In a close-up of a few 
fringes (insets) arrows indicate the field-induced shift. Lower 
panel illustrates the effect of averaging over smooth fluctua- 
tions, Eq. (fT4]h with \a m \ = 2.5s*, which suppress fringe 
contrast without affecting the phase shift. 

shift up in energy by approximately half a period (see 
FigHk). This shift, which is a direct consequence of the 
7r-shift of the reflection phase discussed above, fully de- 
velops in the fields B ~ 0.45*. 

The effect of large-scale potential fluctuations, either 
intrinsic [l6[ or induced by variable distance to gates, can 
be analyzed by averaging the conductance in ([T3]) over 
random offsets in e, taken as a sum of a few harmonics: 



(T) 



Tie — a rn cos(ra0 




(14) 



This simple model describes smooth inhomogeneity with 
the correlation length in excess of p-n interface separation 
L but much shorter than the structure width W. 

The averaging procedure ([T4|) . applied to our numerical 
results, makes the fringes aperiodic and suppresses the 
contrast (see red dashed line in Figj3|). However, the tt 
phase shift induced by magnetic field remains clearly dis- 
cernible even for relatively strong fluctuations (see FigHJ 
where we use J2 m \a m \ = 2.5e* which is quite large com- 
pared to the fringe period of about 0.6e*). 

At even stronger randomness, the FP transmission © 
can be replaced by its phase- averaged value 



(T)e 



tit 2 
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Plugged in ([T3]h it yields magnetoresistance with charac- 
teristic B ~ , identical to that discussed in Q . The re- 
sulting exponential suppression of conductance of course 
holds only in the absence of short-range disorder. 

Conspicuously, the resistance data 11, 12[ feature ape- 
riodic oscillations in gate voltage, observed above the 
point where the sign of carriers beneath the top gate is 
reversed. This is the same region where strong FP fringes 
are expected for an ideal system. The energy scale of the 
oscillations reported in Ref.[ll], converted from gate volt- 
age using 5e/SVtg ~ ^q, is about Se ~ 30 mV, which is 
only a few times larger than the period of 0.86* « 11 mV 
found above (Figs l3|4j) . Could these oscillations, or those 
seen in [l2[ , be the FP fringes contaminated by disorder? 
Comparison to the behavior of the FP contrast in the 
presence of magnetic field, in particular to the tt phase 
shift (FigHJ), may help to clarify this. 

In summary, Fabry-Perot interference in the Klein 
scattering regime is found to be sensitive to the phase of 
the reflection amplitude that exhibits a jump by tt near 
zero incidence angle. This leads to half a period shift of 
interference fringes in the presence of a relatively weak 
magnetic field, a new effect that can help to identify the 
Klein phenomenon in graphene. 
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